Interplay of strongly correlated electrons and localized Ising moments in 

one-dimension 
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We study the ground state properties of the one-dimensional quarter-filled strongly correlated 
electronic chain coupled by J to another chain of antiferromagnetic Ising moments. We focus on 
the case where the large Coulomb interactions localize the charges on every other site. Both the 
electronic spins and the Ising moments interact antiferromagnetically within each chain by J e g and 
J', respectively. Since the number of electrons is half as that of the Ising moments the period of 
magnetic correlation of these two chains are incommensurate. In the presence of J, the frustration 
of J c ff and J' arises, which may lead the system to the intriguing magneto-electric effect. 

PACS numbers: 71.10.Fd, 71.20.Rv, 71.30.+h, 75.50.Dd, 75.47.De 



I. INTRODUCTION 



The electrons coupled to localized spins have been 
a long studied issue from Kondo chains in heavy 
fermionic systemsA, double exchange systems(DEX) in 
manganites^ to 7r-d systems of molecular solids^^. 
Particularly, the negative giant magnetoresistance ef- 
fect in DEX system had been a highlight which pro- 
vided a mechanism of tuning the electronic degrees of 
freedom by the magnetic field^i The author and the 
co-workers tried to figure out another possible mecha- 
nism of the negative magnetoresistance effect by con- 
sidering the Kondo chains including the Coulomb in- 
teraction terms, which we call an extended Kondo lat- 
tice model(EKLM) 8 . The study on EKLM was carried 
out with in mind the experimental findings of the gi- 
ant negative magnetoresistance in the one-dimensional 
(ID) organic solid, TPP[Fe(Pc)(CN) 2 ]2^- This so-called 
phthalocyaninc compound is a quarter-filled 7r-electronic 
system which includes a localized Fe d-spin on the same 
TPP-moleculei^. The 7r-electrons have a charge order, 
which is stabilized by the 7r-d interaction. However, the 
EKLM model turned out to be different from the ex- 
perimental situation mainly in two points; experimen- 
tally, the localized moment has a large Ising anisotropy 
as an artifact of the spin-orbit coupling^, and between 
them, there possibly exists a direct interaction. These 
factors are not considered in the EKLM where the lo- 
calized SU(2) spins do not have the direct interaction 
with each other. We revisit this problem by taking into 
account these points. 

The present paper deals with two chains which have 
electrons and Ising moments, respectively. We focus on 
the strong coupling case where the electrons are localized 
(insulator). Then the system is essentially regarded as 
two spin chains which have antiferromagnetic intra-chain 
correlation with different periodicity. The introduction of 
the coupling between the two chains lead to the "mag- 
netic frustration". The paper is organized as follows; 
Sec. II gives the model and details of the numerical anal- 
yses, and Sec. Ill the ground state phase diagram. Sec. IV 



is devoted to the clarification of the magnetic properties 
of the electrons, and finally the external magnetic field is 
introduced in Sec.V. After all we reach the anticipation 
that the "magnetic frustration" provide a more realistic 
example of the magneto-electric effect. 



II. MODELS AND METHODS 

We consider ID Ising spin chain and electronic chain 
which are connected via Hund coupling — J < at each 
lattice site. Within the chain, the Ising moments interact 
antiferromagnetically by J' > 0, while the electrons have 
strong on-site and inter-site Coulomb interactions, U and 
V, respectively. The Hamiltonian reads, 



TL — Hhubb + Rising + Hj 

Whubb = - 2j(fct o ,Cj .+ h.c. 
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Here, the operators c 
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ij, and Sj denote the annihila- 



tion, number, and the ^-component of the spin at j-th 
site for electrons, respectively. We define the magneti- 
zation of electrons as M c \ — ^ . s| , and the magnetic 
density, m e i = M c \/N, where N is the system size. The 
z-component of the localized Ising moment is represented 
by S? = ±1/2, and (ij) denotes the nearest-neighbor pair 
sites. We focus on the case of quarter-filling of electrons 
where the 4fc^-instability of the charge degrees of free- 
dom is significantly large. 

At J = 0, each of the decoupled chains has well de- 
fined ground state; the Ising chains is a gapped spin 
system with Neel order, which is described as an anti- 
ferromagnetism of 4fc^-periodicity, (tITi ' ' '); where Uf 
is defined by the filling factor of electrons. On the other 
hand, the weak-coupling region of the electronic chain is 
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a Tomonaga-Luttinger liquid (TLL). At quarter-filling, it 
undergoes a transition into a ordered state with a finite 
charge gap 12 ' 13 > 14 when the interactions become as large 
as U > 4t and V •> 2i. In this ordered state at large U 
and V, the electrons localize on every other site. The ef- 
fective antiferromagnetic interaction between the spins of 
these electrons is given as, J® s = t 4 /UV 2 , which leads to 
the correlation characterized as a 2fci?-spin-density-wave 
(SDW). Once J is switched on, there arises a magnetic 
frustration between these two orders. Namely, the anti- 
ferromagnetic configuration of Ising spins favors the elec- 
tronic spins to align ferromagnetically in order to gain 
the energy by J, which is incompatible with the anti- 
ferromagnetic correlation within the electronic chain by 
J c g. Such situation is given schematically in Fig. [Ha). 

Another reference system is the extended Kondo lat- 
tice model (EKLM) in Ref. [8]. This model is realized 
if we replace the interaction of Eq. {T]) between electrons 
and localized spins by that of the SU(2) symmetry, and 
take J' = 0. The role of the quantum fluctuation of the 
localized spins shall be discussed shortly by comparing 
these two models. 

In order to clarify the ground state of this system, we 
perform the density matrix renormalization group anal- 
ysis (DMRG) 1 -. Usually, the Ising anisotropy of spin 
degrees of freedom works as disadvantage in optimiz- 
ing variationally the restricted basis of the finite system. 
This is because the quantum fluctuation which usually 
works to improve the selection of the basis in the local 
update processes is present only in the electron hopping 
term. Therefore, by the straightforward application of 
DMRG to Eq.(JTj), one finds difficulty in selecting the 
optimized configuration of Ising moments. Instead, we 
consider the extended Hubbard chain in the presence of 
periodic field, described by the Hamiltonian, 

n e i = n huhh -J2 w J sZ j- ( 2 ) 

3 

Here, Wj — J(Sj), is the on-site "magnetic field" cre- 
ated by the localized Ising moments. The calculations 
are given as follows; (i) assume several different config- 
urations of Ising moments, (ii) calculate Tt e \ in DMRG 
under the potentials from each of these configurations 
at several system size N, and obtain the energy per site 
in the bulk limit by the finite size scaling analysis, (iii) 
add to (H c \) the interaction energy of Ising moments, 
-Rising = (Rising), and get the lowest energy state as a 
function of J' . As for the Ising spin configuration in (i), 
we consider up to periodicity over sixteen sites, and it 
turns out that the states which have two- or four-fold pe- 
riodicity give the lowest energies. This holds even in the 
non-interacting case, except at J/t < 0.005 where a small 
fcrrimagnetic region is found. We thus mainly focus on 
(TTTT), (TITI), (TUI), and (nit)-configurations, where 
| and | correspond to up and down spin Ising spins, 
respectively. Figures DIa)-[lJd) show the representative 
configuration of Ising moments combined with the unpo- 
larized AF (m e i = 0) and fully polarized F (m e i=0.25) 



electronic states. In the ground state phase diagrams of 
Sec. Ill we consider only the AF- and F-electronic states, 
which allows for the systematic finite size scaling in (ii). 
We confirmed in advance that the partial ferromagnetic 
states with ferrimagnetic Ising moments, e.g. (TTTI), 
do not appear as a ground state in the focus parame- 
ter region. In Sec. IV, we also include the Ising configu- 
rations up to 32-fold periodicity and calculate the m e \- 
dependence by fixing N, in order to examine the effect 
of magnetic field. 

In the finite system with open boundary condition, the 
charges have the largest density at both edge sites. If 
the charge order is dominant the charges tend to align 
in every other sites from both ends, and a calculation 
on even-iV yields a kink structure at the system center, 
e.g. for JV = 8 we find (• o • o o • o •), where • and o 
are the charge rich and poor sites, respectively. In such 
case, the amplitude of the charge density is gradually 
suppressed towards the center site. If we take the odd 
number of sites, this kink disappears, while the electon 
number deviates by one from that of the quarter-filling 16 . 
The results of the finite size scaling of H. c \ up to N = 97 
for (tilt) - and (TTII)-type of Ising moments (which are 
the same configuration in the bulk limit) with both the 
odd- and even-N are presented in Fig. QJb). At N — > oo, 
the energy of all cases coincides within < 10~ 5 £. Thus 
the energy of the bulk limit is safely obtained and the 
treatment (iii) is performed. By introducing finite J', 
the energy density of the (TITI) and (1"f 1"f )-state shifts by 
— J'/4 and J'/4, respectively, whereas that of the (TTII)- 
state does not change. The lowest energy state among 
the calculated candidates are thus obtained as functions 
of U, V, J and J' at fixed t = 1. 



III. PHASE DIAGRAM 

A. Comparison of the ground state of the Ising 
and SU(2) localized moments at J' = 

We first present in Fig. (3Ja) the phase diagram of the 
Ising moments at J' = 0, i.e. when the direct interac- 
tion between Ising moments is absent. This diagram is to 
be compared with the case of the EKLM which has the 
SU(2)-symmetry of localized moments^. At small J, the 
(ftll)-Ising moments couple with the TLL with 2/cp- 
SDW correlation, which undergoes a phase transition 
into the ("|"|"|"|") -Ising moments with fully polarized elec- 
trons. This situation has good correspondence with the 
EKLM; the present (TTII)- and (1"f f|~ )-Ising spin configu- 
rations are interpreted as the paramagnetic and the ferro- 
magnetic state of the SU(2) spins in EKLM, respectively. 
The phase boundary shifts in the Ising case to about four 
times larger value from the SU(2) one. To understand 
this, let us consider the non-interaction case, U = V = 0; 
in the EKLM , the ferromagnetic state at large J is de- 
scribed by the formation of singlet pairs of electrons and 
localized moments, and the rest of the moments (the elec- 
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FIG. 1: (a)-(d) Representative configurations of localized mo- 
ments and electrons in the strong coupling region, (b) En- 
ergy per site as a function of 1/N at U/t = 8, V/t = 4 and 
J/t = 1 with (TUT)- an d (Till)- Ising spin configurations. 
System size scaling is given for four different series; odd N 
with N e = N/2, N e = N/2 + 1, and even TV. 



tron number is half the number of localized moments) 
are fully polarized by the hopping of singlets^. One 
can roughly approximate the energy of the paramagnetic 
and the ferromagnetic states as, e para ~ — 4y2i 



^ and 



eforro ~ ™4t — j, respectively. At J ~ 2t, tne phase tran- 
sition takes place. On the other hand, the energy of the 
corresponding Ising-case is estimated by the modification 
of band structure under the periodic potential of Ising 



2t 



2t 



spins, which yields, e(|f||) ~ 

and e(TTtt) ~ ~4i — j, respectively. The level cross- 
ing occurs approximately at J/t ~ 6, which is consistent 
with the phase diagram. Here, we note that the (Ttll)- 
AF state at small U and V is a band-insulator, because 
the four-fold periodicity leads to the formation of four iso- 
lated bands, where the lowest band is completely filled. 
The above discussions indicates that the lack of quantum 
fluctuation of the Ising moments leads to e(TT-t-L) < e para, 
which originates mainly from the large energy gain of J- 
term in the (tm)-state. Thus the phase boundary is 
pushed towards larger J/t in the present Ising system. 



B. Phase diagram of localized Ising moments 
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FIG. 2: (a) Phase diagram at J' = on the plane oiU/t and 
J/t for V/t = 1,2. Phase boundary of the extended Kondo 
lattice model in Ref. [8] is plotted together for comparison, 
(b) Phase diagrams on the plane of J/t and J'/t for all sets 
of combinations among U/t — 0, 2, 4, 8 and V/t — 1, 2, 3, 4, 5. 
As shown schematically in the inset, the diagrams are sepa- 
rated into four regions, (TTII)-AF, (UTl)-AF, (TITI)-F, and 
(TTTT)-F, which are shown schematically in Figs. l(a)-(d). 
Here, AF and F denote the unpolarized and fully polarized 
electron spins, m e \=Q and 0.25, respectively. 



In the next step, we include the J'-term and find that 
the ground state undergoes a phase transition into an- 
other magnetic state. Figure [3Jb) shows the phase dia- 
gram classified by the configuration of localized moments 
on the plane of J and J' at several fixed values of U and 
V in unit of t. The interaction energies of the Ising mo- 
ments for the configurations (T4T4), (HU), and (TTTT) 
are i?i s i ng /iV = — J'/4, 0, and J'/4, respectively. There- 
fore, by the introduction of J', the ("f J,"f J,)-state which 
gains the energy replaces the others. The boundary of 
the (TITI)- and (TTII)- phases are approximately given 
as J ~ 4 J', which is understood by the comparison of the 
magnetic energy of the (|J,|J,)-AF and (Ttii)-AF states 
in Fig.[TJa); the Ising moments have interaction energies, 
-EWng/A^ = J'/4 and 0, respectively. As for the J-term, 



the straightforward estimation gives, Ej/N = 0(||T|- 
AF) and J/8(HU-AF). However, as we see in Sec. Ill 
(Fig. OHd) at m c \ = 0) , the amplitude of electronic spin 
moment in the (ttjj)-state is (s z ) ~ 0.25, which is sup- 
pressed to about half the expected value. This is pre- 
sumably because the electrons are relatively delocalized 
to the neighboring site (since the neighboring site has the 
same potential from the (ttll)-potential). The resultant 
Ej/N of (||J,J,)-AF is J/16, and after the comparison of 
Rising + Ej in both states, the phase boundary falls on 
J ~ 4J'. 

The phase boundaries are influenced by the electronic 
interactions as well. By comparing the phase diagrams 
we find that U stabilizes the (TTTT)-F state. This is be- 
cause the exclusion of double occupancy due to U fa- 
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FIG. 3: Panels (a) and (b) are the m e i-dependence of energy 
per site, E(m cl )/N, at N = 64 with (Till)- (TfU)-Ising spin 
configuration, respectively, at U/t = 8 and V/t = 4 for various 
choices of J/t. Panels (c) and (d) are the local charge and 
spin densities, {m) and (sf ), which correspond to (a) and (b), 
respectively, at J/t = 0.05. Panels (e) (f) are the spin gap in 
the bulk limit as a function of J/t at U/t — 8 and V/t — 4 
for (Till)- and (TTII)-cases. Panel (g) give the effective spin- 
spin interaction J a a given in Eqs.Q and ((5| as a function of 
J/t at U/t = 8 and V/t = 4. 



vors magnetism. On the other hand, the (TITI)-F state 
(m c \ — 0.25) is stabilized by both U and V. When the 
electronic spins are the fully polarized, J works as poten- 
tials to pin the electrons on every other sites, which favors 
charge order. The (TiTi)-F with large U and V is the typ- 
ical example. Thus U, V, and J cooperatively stabilizes 
the charge ordered state. Such effect of Coulomb inter- 
actions are consistent with what we find in the EKLM 8 . 

The (UU)- and (TTU)-phases at U = 8 given in Fig. 
[UJb) are insulators. We confirmed this by the finite size 
scaling analysis on the charge gap for several choices of 
parameters. The region of U/t > 4 and V/t > 2 is an 
insulator at J=0H. The introduction of J ^ stabilizes 
the insulating phase for both cases, (TIT!) an d (TTID- 
Therefore, our discussion regarding the magnetic prop- 
erties of the insulating state in the next section is safely 
carried out. 



IV. MAGNETIC PROPERTIES 
A. Competing magnetic orders 

The ground state phase diagrams are basically dom- 
inated by the (TIT!) - an d (TTII)- conn g ura tions of 
Ising moments. So far, for each of these cases, the 
spin-unpolarized m c \—0 (AF) and fully spin-polarized 
m e i=0.25 (F) electronic states are examined. In the next 
step, we calculate the m c i-dependence of energy of the 
electronic state by DMRG. Here, we focus on the large U 
and ^-region where the system is a charge ordered insula- 
tor. The nontrivial competition among J, J' and J e g we 
discussed in Sec. II shall be examined explicitly, where J e ff 
is the effective interaction between electronic spins local- 
ized on every other site. Figures [3ja) and [3fh>) show the 
energy per site, E(m e \)/N under the above mentioned 
two different configurations for several choices of J at 
N = 64, U/t = 8 and V/t = 4. In the (TITI)-state, 
the functional form of E(m c \) is flattened as J increases. 
Contrastingly, E(m e i) of the (TTII)-state becomes a more 
rapid increasing function at larger J. 

The charge and electronic spin density for m e i=0 and 
0.25 under (TITI) - an d (TTII)- conn g ura tion are given in 
Figs.[3fc) and^d). The spin densities differ significantly 
between the two figures; at m e i = the (TITI)-state has 
small amplitude of (s z ) compared to that of the (TTII) - 
state. The former has the 4fcp--antiferromagnetic spin 
correlation on the Ising chain which does not fit to that 
of the 2/c^-one along the electronic chain. Therefore, the 
spin density is rather suppressed to cope with J which 
connects the frustrating two different spin correlations. 
In contrast, the latter (TTII)-state has commensurate 
2fci?-antiferromagnetic correlation along both the elec- 
tronic and Ising spin chains. The cooperation of the two 
correlations enhances the amplitude of the spin moment. 
When the electrons become fully polarized at to o i=0.25, 
the magnetic frustration in the (TITI)-state is resolved 
so that the difference between the two configurations be- 
comes almost negligible. 

The spin gap is shown in Figsl3te) and[3jf) as a func- 
tion of J/t at U/t — 8 and V/t — 4, which is obtained af- 
ter the extrapolation to the bulk limit. Both cases starts 
from a gapless state at J—0. Under (TITI)-configuration, 
a small gap opens first and then closes again already at 
extremely small J. As for the (TTII)-case the spin gap 
continues to increase as a function of J. These results 
are consistent with the findings in Figs. 02a) and[3fb). 

As we discussed in Sec. Ill, the charge order is stabi- 
lized by the polarization of the electronic spins in both 
(TITI) an d (TTII)- cases - Actually, the amplitude of (rii) 
increases with increasing m c \. It is interesting to find 
that (rii) does not seem to differ between the (TITI) - an d 
(TTII)- cascs cven at small m e \ where the magnetic prop- 
erties of the two significantly differ. Namely, at quarter- 
filling the charge degrees of freedom is rather decoupled 
to the spin degrees of freedom. This does not hold off- 
quarter-filling which we discuss in Sec. V. 
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FIG. 4: The representative charge ordered state at large U/t 
and V/t for (a) (Ti-Tl)- an d (b)(TT 44)-configurations, where 
the Ising moments work as internal uniform and staggard field 
to the electronic spins, respectively. Panel (c) is the fourth 
order perturbative process by the hopping of electrons. The 
processes themselves are common between (a) and (b) while 
the energy of each state in the processes differ due to different 
configuration of Ising moments. 



B. Strong coupling approach 

In order to understand the nature of such spin degrees 
of freedom in the insulating states we derive the effec- 
tive Hamiltonian of the electronic spins by the perturba- 
tive approach. We start from the strong coupling limit, 
U/t, V/t, U/ J, V/ J — > oo, where the charges localize on 
every other site as shown in Figs. BJa) and HJb) . At 
the first-order perturbative level of J, the 7ii s i ng -term in 
Eq.(p} works as an effective internal magnetic field, H e g, 
on these electronic spins. Then the system is understood 
as the non-interacting SU(2) spins under the (internal) 
magnetic field. The spin-exchange interactions within 
the electronic chain appear in the perturbation processes 
at the fourth order level in terms oit/U, t/V. FigureHJc) 
shows the process which mix the adjacent spins (by two 
lattice spacing). Then, H c \ in Eq.Q is transformed to 
the effective Hamiltonian given as, 



t(l).), (3) 



3—11, (/—integer) 



Here, H c s is a /-dependent internal magnetic field from 
the localized Ising moments. For the ("f J,"f J,)-state we 
have, 



1 



F = — ( 1 

cff ~ U \{V + J/2) 2 + {V - J/2) 2 

, _^ 1 

eff tt M2 I l /l\2 ' 



HcB — 77 



U V 2 - (J/2) 2 
J 



(4) 



where J c z ff > J^ s 



Therefore, the system is an XXZ- 
spin system with Ising anisotropy in the uniform mag- 
netic field. When finite J is introduced the Ising gap 



opens, which, however, is suppressed immediately by the 
magnetic field when J further increases. We show the 
evaluated J e ff's as a function of J/t in Fig. [31(g) from 
Eq.(g]) at U/t = 8, V/t = 4. Actually, the increase of J cff 
is slower than that of -ff e ff- Therefore, even though the 
anisotropy of interaction, J* s /J^, increases as a function 
of J, the uniform magnetic field H c s has larger magni- 
tude and the spin sector is gapless. 

On the other hand, the (ttll)-Ising spin configuration 
yields the following effective parameters, 



J z — 



1 



1 



(V + J/2) 2 J U + J/2 
(1 1,1 



\V 2 (V - J/2) 2 / U - J/2' 
1 1 



V\{V + J){U + J/2) (V — J){U - J/2) 
1 1 



(V+ J/2) 2 {U + J/2) (V — J/2) 2 (U — J/2) / 



(-)'§! J- 



V 2 {V + J/2) 2 ) U + J/2 
~\V 2 + (V- J/2) 2 ) U - J/2 



(5) 



This time we have Jg ff < J^, and the system is inter- 
preted as an XXZ-spin chain with XY-anisotropy (TLL) 
placed under the staggered magnetic field. Figure [U(g) 
shows the increase of J e g as a function of J. Again, the 
internal field overwhelms the effective spin interactions. 
The Neel order is stabilized, which has large spin gap 
compared to the (| || J,)-state. This is simply because 
the period of correlations of Ising and electronic spins 
are commensurate, and they cooperate by the introduc- 
tion of J. The large (s z ) of local moments in the DMRG 
calculation in Fig. [2td) actually supports this scenario. 

We mention that at J = 0, the effective Hamiltoni- 
ans, Eqs.(H]) and are reduced to the simple Heisenberg 
spin Hamiltonian without the magnetic field (i? c ff = 0), 
which has a SU(2) spin-spin interaction, J* s = J^g = 
t 4 /(UV 2 ) ee J c ° ff , which we mentioned in Sec. II. The 
symmetry of J e ff at J ^ is thus modified to a Z(2)-one 
by the localized Ising moments of Z(2) symmetry. 



V. EXTERNAL MAGNETIC FIELD 

A. Phase diagram under external field 

Finally, we introduce the external magnetic field, H , to 
the present system with rather complicated magnetism. 
Figure E^a) shows the phase diagram on the plane of 
J' /t and H/t for several choices of J at N — 64. We 
consider the configuration of local Ising moments up to 
32-site periodicity, all possible degrees of electronic polar- 
ization, m c i, and determine the lowest energy state. The 
phase diagram is classified by the configuration of local- 
ized moments into four parts; at small J' and H ~ 0, the 
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FIG. 5: DMRG results analyzed by the combination with the 
external magnetic field (H) at U/t = 8, V/t = 4 and iV=64. 
Panel (a) is the phase diagram on the plane of J' and H. 
Panel (b) gives the magnetization, m e \ of the (||||)-phase at 
the small H (the corresponding region is indicated by arrows 
above the phase diagram). 



(H ID-phase exists, which is immediately replaced by the 
ferrimagnetic phase by the introduction of H. The sta- 
bilized ferrimagnetic states are the periodic ones, e.g., 
those with four fold, (Till), or eight fold periodicity, 
(TtttttTDi an d finally the full saturation of magnetic 
moments on both chains takes place. The saturation 
field (onset of (TTTT)-F state) of localized moments is 
approximately given as H s ~ 4 J + 2J' . At larger J', the 
(Htl)-phase sustains from if =0 to finite magnetic field. 

Here, again the (||||)-and (fltl)-states significantly 
differ regarding the instability against the external mag- 
netic field. The former state is unstable and disappears 
at finite H. This is because it has no energy gain in 
the Ising part, -Eismg = 0, which allows for the flipping 
of Ising moments into the ferrimagnetic configuration. 
Here, the small spin gap ~ J/2 of the electronic spins 
does not make much advantage. In contrast, the (TITI)- 
phase sustains due to energy gain, Ei s i ng = — NJ'/4, 
which is the largest among all other configurations. In 
this phase, the Ising moment is stable while the mag- 
netization of electrons, m c \ gradually changes with H. 
One can adopt the magnetization curve in Figj5]^b) to 
the (f|f|)-phase ln Fig. [5{a) along the H/t-axis regard- 
less of the value of J'/t. As we saw in Figl^a) the m e \- 
dependence of the energy of this charge ordered state 
was relatively small at finite J. Therefore, the Zeeman 
term of the electronic spins "absolves" the effect of the 
external magnetic field before the localized spins start to 
flip at the J'/t-dependent H = H c (i.e. the (TDD - an d 
ferrimagnetic phase boundary in Fig.[5]Ja)). 



B. Charge gap under the external field 

We discussed in Sec. IV (Fig. [3fc)) that the amplitude 
of charge ordering in the (||||)-state increases by the 
polarization of electronic spins. Therefore, the charge 
gap is also expected to increase with m c \. We calculate 
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FIG. 6: Panel (a) and (c) are "energy gap" at H = (Ao 
without the Zeeman term) and H 7^ (Ah including the 
zeeman term) as a function of m e i and H, respectively, where 
(c) corresponds to the charge gap under the magnetic field. 
The charge gap (Ao at m e i = 0) at H = as a function of 
J is shown together. Panel (b) is the energy of the electronic 
state of (TITI) at N B (= N/2), N e + 1, and N e -1 as a function 
of m c i. Eu, Ev, E.j, and E t are the energies of U, V, J, and 
t-terms in Eq.JTJ. 



the energy gap A = £:(A e -l) + £;(iV e + l)-2£;(iV e ) as a 
function of m e \. For its evaluation at (N e , m e i) with N e — 
N/2, we take either of the magnetic polarization m e i ± 1 
which gives the lower energy in both the electron-doped 
(N e + 1) and hole-doped (N e — 1) states. The result as 
a function of ro e i is shown in Fig. [H^a) for several choices 
of J. As expected, the fully saturated state has larger 
Ao than the unpolarized state. However, at intermediate 
< m e i < 0.25 we find a significant decrease of Ao. This 
decrease is several orders of magnitude larger compared 
to the variation of A induced by J at m c i=0 (see the 
lower panel of Fig. EJa)), and is not attributed to J since 
the decrease is the largest at J = 0. 

Let us consider the origin of the particular m e i- 
dependence of Ao. Figure [B^b) shows the interaction 
and kinetic energies, Ejj, Ey, Ej, and E t , for (iV e )- 
and (A^e ± l)-states, which all together contribute to A . 
What we find is the following; focusing on the energy 
difference between different electron numbers, one sees 
that Eu and Ej do not differ much between (N e )- and 
(N e ± Instates. Whilst, in the (N e + l)-state, Ey and 
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E t have large energy loss at small and large m e i, respec- 
tively, compared to (N e ) and (N e — l)-states. 

Here, E t is an increasing function of m e \, namely the 
paramagnetic state has larger kinetic gain compared to 
the ferromagnetic one. This originates from the Pauli's 
principle as in the usual band picture. At large U and 
V, the doping of electron makes the system less itinerant 
in overall since the electrons try to avoid each other at 
most within the limited spacing. Such tendency is more 
significant for larger m c \ where there is less space for 
electrons to share due to the Pauli's principle. At the 
same time, in the partially polarized case, the electrons 
locally follow the configuration realized at m c \ = to 
gain maximally the kinetic energy. Thus E t at (N e + 1) 
turns out to be a convex downward function. There is 
an energy loss of Ey whenever there is a gain in E t , so 
that Ey in (N e + l)-state behaves contrary to that of 
Et, namely as a convex upward function. Then, these 
convex functions together suppresses the charge gap at 
< TOel < 0.25. 

The convex downward functional form of Ao is the 
most distinct at J = 0, while its both edges have the com- 
parable values, A (m c i = 0) ~ A (m c i = 0.25). When 
J becomes finite, the kinetic energy gain is significantly 
suppressed at m e \ > 0, since the polarized electrons are 
pinned on every other site by the localized moments. 
Thus, E t of the (N e + l)-state at m e i > increases with 
J, particularly at smaller m c \ which was originally more 
itinerant. Then, the convex functional form dissolves and 
Ao is gradually transformed towards a monotonically in- 
creasing function of m e i with increasing J. 

Next, we include the Zeeman term and calculate the 
charge gap against the external magnetic field, which is 
denoted as Ajj- Reflecting the functional form of Ao, 
it appears as a convex-downward function as shown in 
Fig. E^c). Here, for the doped cases we again choose 
either of the m e i ± 1 which gives the lower energy in- 
cluding the Zeeman terms. Notice that we neglect the 
orbital effect under the external field, which we con- 
sider to be small in one-dimension. Particularly at small 
J/t the gap continues to decrease significantly towards 
H/t ~ 0.1 — 0,2. Since the flipping of the Ising mo- 
ments to the ferrimagnetic phase takes place at H c ~ J', 
the upturn of the gap does not appear in overall in the 
(TITD-phaseinFig.^a). 

In this way, the decrease of Ao by m Q \ characteristic 
of the charge ordered state is combined with the stabil- 
ity of the (TITI)Tsing state, and together allow for the 
suppression of charge gap Ah by the magnetic field. 



VI. SUMMARY AND DISCUSSIONS 

In the present paper, we disclosed the intriguing inter- 
play of magnetic and electric properties of the quarter- 
filled strongly correlated electronic chains coupled to the 
Ising moments. Almost regardless of the details of the 
electronic state, the ground state is in overall classified by 



the two different configuration of Ising moments, (TITI) 
and (HID- The latter is an analogue of the paramag- 
netic state of the Kondo lattice model which is stabilized 
by the RKKY-interaction via Kondo coupling, J. When 
the direct interaction between Ising moments, J', is intro- 
duced, this (tt-LJ-)-state is replaced by the (TITI)-state. 

The main focus is the interplay of spin and charge de- 
grees of freedom in the insulating state with large U and 
V coupled to the Ising moments. The electrons are lo- 
calized on every other site, which are interacting anti- 
ferromagnetically via J e g within the chain. The Ising 
moments work as effective internal field |_ff c ff| = J/% 
to these electronic spins. At the same time the Ising 
moments modify the symmetry of the interaction, J ff, 
from the SU(2) at J=0 to the Z(2) ones. Thus, the 
electronic spins coupled to the (T!T!)Tshig moments be- 
have as an Ising XXZ-spin system under uniform mag- 
netic field, and the one coupled to the (TTII)Tsing mo- 
ments is regarded as a XY-spin system under staggerd 
magnetic field. Since H e g is higher in order than the 
Z(2)-modification of interactions, the electronic spins are 
gapless(UU) an d gappcd(TTII), respectively. 

The more simple description is given in the following; 
when the Ising and electronic chains are decoupled, the 
period of magnetic correlation of the electronic chain dif- 
fers by twice from (TITI) while the same as (TTID- hi 
the former case, the magnetic incommensurability of the 
two chains leads to the suppression of the antiferromag- 
netic correlation of electronic spins. In such case, the 
m e i-dependence of energy becomes small, so that the ex- 
ternal field successively flips the electronic spins by the 
Zeeman effect, while the (TITI)Tsing order sustains due 
to the large energy gain by J'. Along with the polariza- 
tion of the electronic spins, the charge gap decreases sig- 
nificantly by H . Contrastingly, the (TTII)-configuration 
is stabilized solely by the magnetic commensurability of 
two chains (and not by J'), thus is easily destroyed by 
the external magnetic field. We conclude that the (TITI) - 
order of Ising spins appear as a consequence of the direct 
interaction between Ising moments J', where a frustra- 
tion of competing magnetic orders is embedded. Such 
frustration yields a strong-corrclation-driven interplay of 
spin and charge degrees of freedom. 

The similar picture shall also be found in the EKLM. 
The paramagnetic state of EKLM corresponds to (TTII) 
in the present study. The (TTII)-ordcr is weakened in the 
EKLM by the quantum fluctuation and transformed into 
the antiferromagnetic correlation which cooperates with 
Jeff via RKKY interactional. When J becomes large, the 
SU(2) spins form singlets with the electrons, which prop- 
agate and stabilize the ferromagnetism- 7 . Therefore, if 
one includes the direct interaction, J', between SU(2) 
spins, it favors antiferromagnetism and competes with 
the J- induced ferromagnctism. Thus, a similar mag- 
netic frustration may appear. However, such physical 
picture shall be rather blurred by the quantum fluctua- 
tion(SU(2)). The present Ising spin system has a more 
serious frustration effect, which may lead the electronic 
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system to a sensitive response against the magnetic field. 

Finally, let us examine the relevance of these re- 
sults with the experiments on the TPP[M(Pc)(CN) 2 ] 2 , 
M=Fe,Coi£. The Co-salt is a pure electronic chain 
(J = in the present model), which is a good reference 
system to analyze the effect of localized moments. Both 
salts have semi-conducting temperature dependence of 
resistivity. The Fe-salt shows a large negative magne- 
toresistance (MR) which amounts to p(H)/p(0) ~ 10~ 2 , 
where p{H) is the resistivity under the magnetic field 
H. The experimental findings are summarized as fol- 
lows; (1) the activation energy derived from the resis- 
tivity data are A a ~ 10 _3 eV and ~ 10 _2 eV for Co-and 
Fe-salts, respectively^. (2) The localized moment has 
anisotropic g-values, g± ~ 3.6 and gu ~ 0.5 — 1 which are 
roughly perpendicular and parallel to the molecular-axis, 
respectively. This anisotropy is explained in terms of the 
spin-orbit coupling and the resultant magnetic moment 
is S = 1/2 11 . (3) The magnetic susceptibility shows a 
large anisotropy, X±/X\\ > 5, and \\\ gives similar values 
with the Co-one^. (4) Residual magnetization is observed 
at T < 12K, which is attributed to the ferrimagnetism 
of 7r-electrons by the torque experiment^. (5) The MR 
is not scaled by the magnetization and shows large In- 
dependence^. (6) Phase transition is absent (in con- 
trast to DEX), namely the scenario of the competition 
of two different orders are not applicable. (7) The MR 
is relevant when Fe-ion is partially replaced by Co-ion, 
e.g. even when the degree of replacement is as large as 
Feo.07Coo.93, the MR amounts to p(H)/p(0) ~ 0.522, (8) 
the magnetization of electrons gradually increases with 
H, and that of the Fe-moments starts at H c ~ 15T 21 . 
(9) the ground state of the electronic chain is weakly 
charge ordered which is observed by the NQR studyi^. 

The factor (2) and (9) are taken into account in the 
present model. The parameter values estimated from the 



extended Hiickel calculation, ab-initio calculation and by 
the reflectance spectra give t ~0.1eV, J ~ t/10, and 
J' ~ J/3, where J and J' are Hund and exchange cou- 
pling constants, respectively. If we put these parameters 
onto the phase diagram of Fig.[2jb), it is located within 
the (tit!) -phase. The magnetic properties of electrons 
behave quite sensitive to U and V; when U = 8 and 
V — 4, we see that the system is in the vicinity of the 
phase boundary between the AF- and F-electrons. As we 
saw in Fig. [3]Ja) the m c i-dependence of energy is small. 
Then, the ferrimagnetism of 7r-electrons is possible, which 
is consistent with (4). The external field gradually mag- 
netizes the 7r-electrons, and the onset of flipping of local- 
ized moments, H c ~ O.Olt in (8) is consistent with the 
phase digram in Fig. [5Ja) (at J'/t~0.1/3). In this way, 
we reach the picture that electrons show fragile magnetic 
properties under the stable antiferromagnetism of local- 
ized Ising moment supported by J'. 

However, there still remains some issues to be clarified; 
in our model, the charge gap decreases with H, while 
it does decrease even at J = 0. Therefore, J does not 
seem to favor the suppression of charge order by the mag- 
netic field. While, the Co-salt shows a small but positive 
MRji2, which is incompatible with our results. Also the 
present model which focuses on the particular density of 
d-spins cannot cope with the issue (7). Therefore, in or- 
der to clarify fully the origin of magnetoresistance, a more 
systematic experimental data (regarding the Fe-Co-ratio 
or the effect of the dimensionality of the system) , as well 
as the theoretical calculations directly on the transport 
properties shall be required. 
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